Order parameter sign-reversal near s±-superconductor surface 
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The superconducting order parameter and LDOS spectra near an impenetrable surface are studied 
on the basis of selfconsistent calculations for a two band superconductor with nodeless extended s- 
wave order parameter symmetry, as possibly realized in Fe-based high-temperature superconductors. 
It is found that for a wide range of parameters the spatial behavior of the order parameter at 
a surface is not reduced to a trivial suppression. If the interband scattering at a surface is of 
the order of the intraband one or dominates it, it can be energetically favorable to change the 
symmetry of the superconducting state near the surface from s± to conventional s-wave. The range 
of existing this surface conventional superconductivity is very sensitive to the relative values of 
interband and intraband pairing potentials. It is shown that the LDOS spectra near the surface can 
qualitatively differ upon calculating with and without taking into account the selfconsistency of the 
order parameter. 

PACS numbers: 74.45.+C, 74.70.Xa 



o 

o 
1 

Sh 1 

3 : 

a ■ 



-a 
c 

o 
o 



> 

CO 
CN 

co 

o 
o 



X 



The discovery of a new family of iron-based high- 
temperature superconductors with distinct multi-orbital 
band structure* — has renewed interest to the problem of 
multi-band superconductivity, firstly discussed fifty years 
ago^£. It was proposed theoretically 6 -! 7 - that the Fe-based 
superconductors represent the first example of multigap 
superconductivity with a phase difference between the su- 
perconducting condensates belonging to different bands. 
This state was discussed previously^ 9 , but not yet ob- 
served in nature. In the most simple case there is the 
phase difference 7r between the superconducting conden- 
sates arising on the hole Fermi surfaces around T point 
and the electron Fermi surfaces around M point. This so- 
called s± (or extended s-wave) state has been favored by 
a variety of models within random phase approximation 
(RPA p 10 ' and renormalization group techniques^—. 
Currently the s±-state is viewed to be the most plausi- 
ble candidate for the role of the superconducting order 
parameter in these compounds. 

Surface and interface phenomena in s±- 
superconductors have attracted considerable recent 
attention. The formation of bound states at a free sur- 
face of an s±-superconductor— ~— , at an S±/Ni2r— , an 
N/S/S± junctional and at Josephson junctions including 
s±-superconductor s 22 ' 24 was investigated theoretically. 
In particular, the finite energy subgap bound states 
(depending on the interface parameters) were found 
and their influence on the conductance spectra and 
Josephson current was investigated. 

However almost all these calculations (except for a 
few numerical results^) assume non-selfconsistent super- 
conducting order parameter (OP). In the present paper 
we focus on the study of the OP at a surface of s±- 
superconductor. We have found that for a wide range 
of parameters the spatial behavior of the OP at a sur- 
face can not be reduced to a trivial suppression. If the 
interband scattering at a surface R12 is of the order of 
the intraband one Rq or dominates it, it can be energeti- 
cally favorable to change the symmetry of the supercon- 



ducting state near the surface from s± to conventional 
s-wave. The range of existing this surface conventional 
superconductivity is very sensitive to the relative val- 
ues of interband and intraband pairing potentials. We 
demonstrate that the selfconsistent OP behavior affects 
the surface LDOS profiles, and, consequently, should be 
taking into account when interpreting experimental re- 
sults. It is worth to note here that, while there is a wide 
parameter range of existing complex OP at the surface 
region—, in this paper we only discuss the case when the 
surface OP is of conventional s-wave type. 

We consider an impenetrable surface of a clean two- 
band superconductor. The OP is assumed to be of s±- 
symmetry in the bulk of the superconductor, that is 
the phase difference between the OP's in the two bands 
(called 1 and 2) is n. It is supposed that an incoming 
quasiparticle from band 1,2 can be scattered by the sur- 
face as into the same band (intraband scattering), so as 
into the other band (interband scattering). 

We make use of the quasiclassical theory of supercon- 
ductivity, where all the relevant physical information is 
contained in the quasiclassical Green function gi(e,pt, x) 
for a given quasiparticle trajectory. Here e is the quasi- 
particle energy measured from the chemical potential, pf 
is the momentum on the Fermi surface (that can have 
several branches), corresponding to the considered tra- 
jectory, x is the spatial coordinate along the normal to 
the surface and i = 1, 2 is the band index. Quasiclassical 
Green function is a 2 x 2 matrix in particle-hole space, 
that is denoted by the symbol The equation of mo- 
tion for gi(e,pf,x) is the Eilenberger equation subject to 
the normalization conditio n 26 : 27 . For superconductivity 
of s±-type, when the pairing of electrons from different 
bands is absent, the Eilenberger equations corresponding 
to the bands 1 and 2 are independent. The trajectories 
belonging to the different bands can only be entangled 
by the surface, which enters the quasiclassical theory in 
the form of effective boundary conditions connecting the 
incident and outgoing trajectories. 
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However, owing to the normalization condition for the 
quasiclassical propagator, the boundary conditions for 
the quasiclassical Green functions are formulated as non- 
linear equations 2 ^— . Furthermore, they contain unphys- 
ical, spurious solutions, so their practical use is lim- 
ited. For this reason in the present work we make use 
of the quasiclassical formalism in term of so-called Ric- 
cati amplitude o 31 ' 32 , that allows an explicit formulation 
of boundary conditions 3 ^—. The retarded Green func- 
tion gi(e,pf,x), which is enough for a complete descrip- 
tion of an equilibrium system, can be parametrized via 
two Riccati amplitudes (coherence functions) ji(e,pf,x) 
and ji(e,pf,x) (in the present paper we follow the no- 
tations of Refs . 32 ' 36 ). The coherence functions obey the 
Riccati-type transport equations. In the considered here 
case of two-band clean s±-superconductor the equations 
for the two bands are independent and read as follows 



A7i + 2£7* = -A* 7l 2 
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Here Vi x is the normal to the surface Fermi velocity com- 
ponent for the quasiparticle belonging to band i. Aj 
stands for the OP in the i-th band, which should be found 
self-consistently. 

Let us suppose that the surface is located at x = 
and the superconductor occupies the halfspace x > 0. 
For the sake of simplicity we assume that the surface 
is atomically clean and, consequently, conserves paral- 
lel momentum component. Then there are four quasi- 
particle trajectories, which are involved in each surface 
scattering event. These are two incoming trajectories be- 
longing to the bands 1,2 (with v ix < 0) and two outgoing 
ones (with Vi X > 0). It can be show n 32 ' 36 that the coher- 
ence function 7i(e,p/,x), corresponding to the incoming 
trajectory can be unambiguously calculated making use 
of Eq. ([1} up to the surface starting from its asymptotic 
value in the bulk 
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where A\ is the bulk value of the OP in the appropriate 
band, d > is an infinitesimal. As for the coherence 
function ji(s,pf,x), it is determined unambiguously by 
the asymptotic conditions for the outgoing trajectories 
and can be obtained according to Eqs. ((TJ),©- 

Otherwise, the coherence functions ji(e,pf,x) for the 
outgoing trajectories and, correspondingly, ii{e,pf,x) 
for the incoming ones should be calculated from Eq. ^ 
supplemented by the boundary conditions at the surface 
and Eq. ((2]) . The surface is described by the normal state 
scattering matrix for particle-like excitations, denoted by 
S and for hole-like excitations, denoted by S, that con- 
nect outgoing with incoming quasiparticles. The scatter- 
ing matrix S have elements Sk iPj , which connect outgo- 
ing quasiparticles from band i with momentum hi to the 



incoming ones belonging to band j with momentum pj . 
Here and below all the momenta corresponding to the 
incoming trajectories are denoted by letter p and all the 
momenta for the outgoing quasiparticles are denoted by 
fe. For the model we consider S is a 2 x 2-matrix (for the 
particular value of the momentum parallel to the surface) 
in the trajectory space. It obeys the unitary condition 
SS^ = 1 and without loss of generality can be parame- 
terized by three quantities -R12, 6 and a as follows 
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where i?o and R\2 are coefficients of intraband and inter- 
band reflection, respectively. They obey the constraint 
R0 + R12 = 1- The phase factors a = ±1 and O appear to 
be unimportant for further consideration. While in gen- 
eral the scattering matrix elements are functions of the 
momentum parallel to the surface pn, we disregard this 
dependence in order to simplify the analysis. The scat- 
tering matrix S for hole-like excitations are connected 
to S by the relation S(p\\) = S tr (— pn). In the absence 
of spin-orbit interaction the S'-matrix elements are only 
functions of that is in the case we consider S = S. 

From the general boundary conditions 3 - 6 -, which arc 
also valid for a multiband system, one can obtain the 
explicit values of the coherence functions 7,(£, fe, x = 0) 
and 7i(e, p, x = 0) via the scattering matrix elements and 
the values of the coherence functions ji(e,p,x = 0) and 
7i(e, fe, x = 0) at the surface. They read as follows 



Tife = -Ro7ip + -Ri272p- 



flo-Ri272fc(7i P ~ 72 P ) 2 

1 + 72fc (-Rl27lp + -Ro72 P ) 



lip = RQllk+Rl2l2k-— 7^—3 , p ~ s 

1 + 72p (itl27lfc + Rol2k) 



, (5) 



(6) 



Here the arguments (e, x = 0) of all the coherence func- 
tions are omitted for brevity, j ip = 7, (p) and 7j p = 7, (p) 
and the analogous notations are used for 7i(fe) and ji(k). 
Quantities 7.^ and 7^, entering Eqs. ([5]) and (J6j) are to be 
calculated from Eqs. (JTJ and ([2]) supplemented by the ap- 
propriate asymptotic condition. The coherence functions 
72fc and 72 P are obtained by the interchanging 1 f-> 2 in 
all the coherence function band indices at the right-hand 
side of Eqs. §5§ and ©, respectively. 

Now, substituting the coherence functions into the self- 
consistency equation 



A t (x) = -T^Xi 



-2Z7T7 



3Pf 



1 
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(7) 



we iterate system ([I])-©, ©-(0 until it converges. In 
Eq. Xu < is the dimensionless pairing potential 
for band i and A12 = A21 is the dimensionless interband 
pair-scattering potential. We choose A12 > 0, which 
stabilizes s± OP in the bulk. The Matsubara frequen- 
cies e n enter the coherence functions via the substitution 
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e + iS — > ie n . (■■■}pf means the anomalous Green function 
averaged over the entire Fermi surface, that is Pf incor- 
porates as the incoming trajectories p, so as the outgoing 
ones k. For concreteness we suppose the Fermi surface 
to be cylindrical for the each band. However, our results 
do not qualitatively sensitive to this assumption. 
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FIG. 1. Spatial profiles of the OP corresponding to 
band 1 (panel(a)) and band 2 (panel(b)) near the surface 
(x — 0) for different values of interband pair scatter- 
ing. R 12 = 1, Ajj = -0.25A?, temperature T = 0.3A?. 
The values of intraband pairing potentials \a are ad- 
justed to keep A\ and A 2 unchanged upon varying 
A12. The particular values of coupling constants are 
the following: (Ai 2 , An, A 22 ) = (0.06,-0.2607,-0.0327) 
for black solid curve, (0.03,-0.2695,-0.1353) for 
dashed curve, (0.01, -0.2753, -0.2036) for dotted curve, 
(0.004,-0.2771,-0.2241) for dashed-dotted curve and 
(0.002, -0.2777, -0.2310) for gray solid curve. The super- 
conducting coherence length £ = Vi// A\ . 

The spatial profiles of the OP calculated according to 
the described above technique, are represented in Figs.QJ 
H We assume that in the bulk |Aj| > |A£|. Panels (a) 
of Figs. HJ3] demonstrate the spatial OP profiles for band 
1, while panels (b) correspond to band 2. 

Fig. []] shows the dependence of the effect on the in- 
terband pair-scattering value A12. While the larger OP 
(in band 1) is simply suppressed near the surface and 
the magnitude of the suppression is only slightly sen- 
sitive to A12 (at least, in the range we consider), the 
smaller OP (in band 2) reverses its sign near the surface. 
Thus, there is a surface region, which size is comparable 
to the superconducting coherence length £, where s±- 
superconductivity is superseded by the conventional one. 
The reason for this OP sign reversal is the interband 
surface scattering R 12 , because it is energetically more 
favorable to minimize the OP gradient term along the 
quasiparticle trajectory. It is worth to note here some- 
what similar effect of the OP sign reversal for the smaller 
gap due to magnetic impurities^. 

Fig. Q] demonstrates that the weaker the interband pair 
scattering potential the wider the region, where the con- 
ventional s-wave superconductivy exists. It can be qual- 
itatively understood on the basis of the self-consistency 
equation (J7J: the phase of the smaller OP is determined 
to a great extent by the phase of the dominant OP via 
the term A12. The weaker this connection the less en- 
ergy cost to reverse the phase. As it is seen from the 
caption to Fig. [TJ in the particular calculations we take 
I A12 ] << I An I - Such a choice of parameters is consistent 



with the experimental estimates of the coupling constants 
for FeSe2I. 
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FIG. 2. Spatial profiles of the OP near the surface in band 
1 (panel (a)) and in band 2 (panel (b)) for different values of 
interband surface scattering R 12 . a\ = -0.25A?, T = 0.3A?, 
(A 12 , An, A 22 ) = (0.004, -0.2771, -0.2241). R 12 = 1 for black 
solid curve, 0.8 for dashed curve, 0.5 for dotted curve and 0.2 
for dashed-dotted curve. 

The curves represented in Fig. [T] are calculated un- 
der the assumption of purely interband surface scattering 
i?i2 = 1, when the OP sign reversal is strongest. In or- 
der to investigate the effect in the more realistic situation 
one needs to take into account intraband scattering Rq. 
The corresponding results are demonstrated in Fig. O It 
is seen that the region of reversed OP existence shrinks 
with increasing of Rq. However the effect remains to be 
pronounced even if the portion of intraband scattering 
exceeds 50%. Thus, we believe that the self-consistent 
OP treatment is essential for polycrystalline samples, for 
example, upon analysing spectroscopic data. 
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FIG. 3. Spatial profiles of the OP near the surface in band 1 
(panel (a)) and in band 2 (panel (b)) for different values of the 
ratio AS/ Aj. T = 0.3A?, A12 = 0.005, R12 = 1. A 2 /A? = -1 
(An = A 22 = 0.2732) for gray solid curve, AS/ A? = -0.5 
(An = 0.2754, A 22 = 0.23 93) for black solid curve, A 2 /A? = 
-0.25 (An = 0.2768, A 22 = 0.22 07) for dashed curve and 
A|/A? = -0.125 (An = 0.2775, A 22 = 0.2011) for dotted 
curve. 

The considerable region of conventional s-wave surface 
superconductivity can only occur if the bulk OP's in the 
two bands essentially differ in magnitude. In case if |A|| 
approaches to |Aj|, the OP phase reversal region shrinks 
and, finally, the surface OP behavior reduces to the triv- 
ial suppression for the two bands if \A\\ = |A||. Oth- 
erwise, upon decreasing the ratio |A||/|At the surface 
region gets wider until its width saturates at some value 
of the ratio (~ 1/4 for the considered case). This is illus- 
trated in Fig. [3] 

Now we discuss how the described above OP spatial 
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FIG. 4. LDOS as a function of quasiparticle energy calcu- 
lated at x = 0. The energy is measured in units of A*. 
Left column represents LDOS for band 1, while right col- 
umn corresponds to band 2. Upper row demonstrates the 
results of selfconsistent calculations and lower one shows 
the LDOS assuming non-selfconsistent OP. (A12, An, A22) = 
(0.004, -0.2771, -0.2241), A2 = -0.25A5, T = 0.3A?, R 12 = 
0.5. 
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FIG. 5. The same as in Fig. [4] but corresponding to R12 = 1. 

behavior affects the local density of states (LDOS) near 
the surface. LDOS pi corresponding to band i is calcu- 
lated via the coherence functions as follows 

p 4 (*, £ )=Re/i-^4^\ . (8) 
\1 + 7<P/7<P/ / Pf 

Figs.|4]and[5]represent the LDOS at the surface {x = 0) 
as a function of quasiparticle energy. Left and right 
columns of the Figures demonstrate the LDOS for bands 
1 and 2 separately. Upper row of each Figure shows 
LDOS plots, calculated taking into account the selfcon- 
sistent OP behavior. It should be compared to the lower 
row, where LDOS is plotted for the same parameters, 
but for the non-selfconsistent OP equal to its bulk value. 
The results represented in Fig.0]correspond to R12 = 0.5, 



while Fig. [5] illustrates the case of purely interband scat- 
tering R12 = 1. 

As it was already discussed in the literature, if R12 7^ 
there are surface bound states in the system, which man- 
ifest themselves as well-pronounced peaks in the LDOS. 
At i?i2 — > 1 the bound state energies tend to zero, what 
can be clearly seen from the corresponding LDOS plots. 
For this case the LDOS is dominated by very strong zero- 
energy peak and the differences between selfconsistent 
and non-selfconsistent plots are not qualitative. It can 
be only noted that the features corresponding to gap 
edges (clearly seen at least in panel (c) of Fig. [5]) are 
washed out under selfconsistent calculation. At the same 
time for the intermediate value of interband scattering 
there are qualitative differencies between selfconsistent 
and non-selfconsistent results (see Fig. 0} . They can be 
summarize as follows: (i) while in the non-selfconsistent 
picture the bound state peaks are divided by the clearly 
defined gap, selfconsistency results in transforming this 
inner gap into "V" -shaped behavior, which is known to 
be more typical for the superconductors with OP nodes 
at the Fermi surface. We believe that this observation 
can be essential for interpreting experimental data, (ii) 
additional features (small peaks) appear in the subgap re- 
gion upon taking into account selfconsistency. It is worth 
to note here that, in contrast to the interface OP behav- 
ior, the shape of LDOS profiles can be quite sensitive to 
the details of the microscopic model describing the in- 
terface, in particular to the concrete dependence of the 
scattering matrix elements on pn. However, if the par- 
ticular microscopic scattering matrix model leads to the 
existence of an inner gap in the LDOS, it is inevitably 
transformed into "V" -shaped behavior under selfconsis- 
tent calculation, as it is demonstrated above. This fact 
is a consequence of spatial line of OP nodes appearing if 
the OP sign reversal takes place at the surface. 

In summary, for a two band s±-superconductor we 
have theoretically investigated the behavior of the OP 
at a specular reflecting surface. It is found that if the 
interband surface scattering is of the order of the intra- 
band one or dominates it, the OP belonging to the band 
with smaller OP absolute value in the bulk reverses its 
sign at a surface region, thus giving rise to conventional 
s-wave surface superconductivity. This region of reversal 
sign OP is maximal for purely interband surface scat- 
tering and shrinks upon its diminishing. The effect is 
quite sensitive to the ratio between intraband and in- 
terband superconducting coupling constants and is more 
pronounced for smaller interband coupling constant. It 
is also shown that if the OP sign reversal takes place, it 
results in qualitative changes in the LDOS spectra near 
the surface. 
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